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Abstract 

The equation of state and composition of matter are calculated for conditions typical 
for pro-collapse and early collapse stages in core collapse supernovae. The composition 
is evaluated under the assumption of nuclear statistical equilibrium, when the matter is 
considered as an 'almost' ideal gas with corrections due to thermal excitations of nuclei, to 
free nucleon degeneracy, and to Coulomb and surface energy corrections. The account of 
these corrections allows us to obtain the composition for densities a bit below the nuclear 
matter density. Through comparisons with the equation of state (EOS) developed by Shen 
et al. we examine the approximation of one representative nucleus used in most of recent 
supernova EOS's. We find that widely distributed compositions in the nuclear chart are 
different, depending on the mass formula, while the thcrmodynamical quantities are quite 
close to those in the Shen's EOS. 

keywords: Nuclear reactions, nucleosynthesis, abundances - Stars: neutron - super- 
novae: general 



1 Introduction 



A proper description of the equation of state (EOS) for subnuclear and supernuclear den- 
sities is of v ital importance for current studies on the ex plosion mechanism of core-collapse 
supernovae (|Bethelll99d ISuzuki|[l99i Ijanka et alll2007l) . The success of the prompt shock 
propagation in the supernova mechanism depends on the size of the homologously contract- 
ing inner part of a collapsing core. The larger is Yc (the number of electrons per nucleon) 
at the time of bounce, the smaller part of matter will be dissociated and the stronger will 
be the prompt shock wave. Therefore, all factors which influence Y^. must be thoroughly 
studied. 

One factor is the mass fraction of the free protons, Xp, in the dense matter of supernova 
core. Since the rate of electron capture on a free proton is much higher than that on a 
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nucleus, free protons play an important role in establishing the value of Ye at the neutrino 
trapping. The abundance of free protons depends sensitively on the nuclear models of dense 
matter. The variati ons of arnount to more than order of magnitude in the studies of su- 



pernova EOS so f ar (ICooperstein "l985: C ooperstein fc Baron 1990:' Hillebrandt fc Woll l985t 
Hillebranddll991 : iLattimer et al.iil985, : ,Lattimer fc SwestVMl991i:iShen et al.lll998allblV and can 



affect the i nitial strength of shock wave |Bruennl 1989a : Swestv et al. 1994 : Sumivoshi et al . 
2004l2005l) . As shown bv lBruennI (|l989jbl ). the free proton fraction may be largely different 
depending on the nuclear interaction and its model. Another factor is the composition of 
various nuclei, which exist under the nuclear statistical equilibrium. The total rate of elec- 
tron cap tures on nuclei c an dominate that on free protons, if the number of the free protons 
is small ( Hix et al. and nuclei may affect the dynamics in this way as well. Therefore, 

it is crucial to evaluate the composition of dense matter in a precise manner. 

Another astrophysical problem, where an accurate chemical composition is important, 
is the nucleosynthesis of heavy elements. The studies of various m odels of the explosive 
nucleosynthesis as well as the rapid neutron capture process Csee. e.g. iPtitsvn fc Chechetkin 
1982t IWooslev fc HoffmaJll992t iNadvozhin et al.l Il998t lArnould et al.l 120071 and references 
therein) require knowledge of the chemical composition close to the nuclear statistical equi- 
librium (NSE) as an initial condition for nuclear reaction network calculations. Neutron star 
mergers and collapsar models for gamma-ray bursts may contain similar conditions as well. 

The thermodynamic properties of h ot, dense matter ha v e been investigated in var- 
ious approaches: Saha-like equati o ns ( Mazurek et al. _ 19791: IeI Eid fc Hillebrandtl Il980l: 



Murph-y 



2proacn( 
\lll980l: 



Ishizuka et a l. 2003; Nadvoz hin fc YudinI 12004) . Hartree-Fock appr oaches 



press ible liquid-drop model ( Bavm et al.l Il97ir Lattimer et al. Il985t Lattimer fc Swestv 



(jBonche fc Vautherinlll98 ir iWolfil983;.HillebraiKlt fc Woljll9'8a iLassaut et al.lll987h 



com- 



1991 ) , a nd the relativistic m ean field theory ( Sutaria et al.l 19991 ) with local de nsity approx- 



imatio n (jShen et al.lll998allbl ). The review of those approaches, be ginning from | Bethe et al 



(Il970l). can be found in the context of cold neutron star matter in iRiister et al] (|2006f) and 



Haensel et al.l (|2007n . In most of the sets of EOS used for recent supernova simulations, the 
approximation of o ne species of nuclei ( i n addition to neutrons, protons and alpha particles ) 
has been adopted (|Hillebrandt fc Wold [19851 : ILattimer fc Swestv! [l99ll IShen et al.lll998bh . 
Howev er, the advance in recent studies of electron capture rates on nuclei for a wide mass 
range ( Langanke fc Martfnez-Pinedol 2003 ) necessitates the evaluation of multi-composition 
of nuclei to determine the total electron capture rate. Although the treatment with one nu- 
clear species sho uld be a good approximatio n to derive overall thermodynamical properties 
of dense matter ( Burrows fc Lattimer 1984 ). it is necessary to take into account the multi- 
composition at high densities ~10^^ g/cm^ near the neutrino trapping regime, where nuclear 
interaction and multi-compositional treatment may become more influential to determine 
the equation of state. The description of the EOS with the multi-c omposition is, o f course, 
crucial to predict the abundance of nuclei and free protons, see e.g. iMurohvl (Il98fll ). 

In spite of an extensive use of the Saha equation in the literature, not all relevant 
physics was ii icluded i n the p ublished work within this approach from the very beginning. 
For example, Murphv (Il980l) has not taken into account the effects of non-ideal nucleon 
gas. IeI Eid fc Hillebrandtl (|l980h have consistently ta k en int o account finite temperature 
effects in nucleon interactio ns following El Eid fc lilj (Il977[). but they have omitted the 
Coulomb corrections. Lat er Hillebrandt fc Miilled ( 198l[ ). iBravo fc Garcia-Sen 3 l|l999h and 



Nadvozhin fc YudinI (|2005l ) included them in different app roximations and solved the implicit 
set of Saha equations. Moreover, Hillebrandt fc Miillerl (|l981. ') have tak en into account an 



'excluded- volume' effect due to the finite size of nuclei. It was noted bv lHillebrandtl ()l99lh 
that the number of nuclear species included in the Saha treatment strongly affects the mass 
fraction of free protons, which may change the effective adiabatic index through electron 
captures. 

In this paper we describe our code suitable for studying the properties of dense matter 
relevant to core-collapsing supernova conditions, and we put an emphasis on the composi- 
tion of nuclei. Our code is originally rather old: we have compared our NSE results with 
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the kinetic approach aheady in the paper ( Panov et al.ll200"ll ). but the code has not been 
desc ribed in detail. Ano ther motivation for our w ork is the availabihty of new ma ss formulae 
like (iKoura et al.ll2005h and (|M511er et al.lll995h and detailed EOS tables like (|Shen et al 



1998a) . Our aim is to develop and to test a practical and reliable tool for predictions of NSE 
composition at subnuclear densities. 



Followi ng most closely t he co nventional Saha approach (jClifford fc Tavleij|l965f ) and the 
method bv lMazurek et all (|l979l ). we extend this approach in the following points: 



We take into account, at a certain level of approximation, the influence of free nucleon 
gas on the surface and Coulomb energies of nuclei. We retain some terms which were 
omitted bv lMazurek et al. I (ll979l) . However, we do not take into account the 'excluded- 
volume' effect, not pretending to reach very high densities. 

W e have an opt i on to include various results for nuclear p artition functions like those 
of Fowler et al. I (Il978l) and newer p artition functions by lEngelbrecht fc Engelbrech^ 
(|l99ll ) and lEngelbrecht eTal] (|l990li . 



Our network is appreciably larger than the one used previously. The atomic mass table 
is updated using recent theoretical compilations of at omic masses. It co vers ^^20000 
nuclides (Koura 2007) for the KTU Y mass formula (iKoura et all |2005| ) and ~9000 
nuclides for the FRDM mass formula (jMoller et al.lll995l ) as an extra option. 



To examine the basic properties of the equation of state for supernovae in the current 
formulation, we calculate the properties of dense matter covering a wide range of density, 
p, electron fraction, Y^, and temperature, T . We report here the results at the equilibrium 
for the nuclear and electromagnetic processes, i.e. the NSE, for fixed values of without 
impo sing beta equilib r ium. The b eta-equilibrium is easil y incorporated in the Saha approach 



Mazurek et al] (|l979t ) and IeI Eid fc Hillebrandtl (fl98Q ). and our code has such an 



option. A brief discussion of account of the beta equilibrium is also given below. We show 
the global features of the dense matter (thermodynamical quantities and compositions) for 
the supernova environment. In order to assess the dependence of the nuclear mass, we make 
the comparisons among several choices. We also examine the similarity and diff erence from 



the sin gle nuclear species treatment by comparisons with the Shcn EOS table IjShen et al 
19983). 

We explain the formulation of nuclear statistical equilibrium in §2.11 with the detailed 
description of spin of nuclei in §2.21 We describe briefly the atomic mass data used for the 
calculations in §2.31 the Shen equation of state in §2.41 and the beta equilibrium in §2.51 
We show the properties of dense matter in the current framework in fJ3) The summary and 
discussions will be given in Sj4l 



2 Thermodynamics of interacting nuclides 

2.1 Inclusion of nuclear and Coulomb contributions 

The basis is the equation of nuclear statistical equilibrium (NSE) with respect to strong and 
electromagnetic processes for the chemical potentials of nuclei: 

Ml = ■^iMp + {^i - Zi)p,n ■ (1) 

We denote i — 1 = p for proton and i = 2 = n for neutron. We use the index i mostly 
for nuclei (Ai, Zi) with A > 1, but we include nucleons (i = 1, 2) in the summation over all 
species. To find the correct expressions for fXi we suppose that the free energy of the set of 
{Ni} nuclides is 

^^„uc = ^iV,$,(r,{nJ) . (2) 
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Here $i is the free energy per nucleus i and it may depend not only on the concentration, 



(3) 



of the i-th nucleus but also on concentrations of some other nuclei rik- Then we have by 
definition 



dN, 



T,V 



E 



nk 



dn,. 



Thus we find for the pressure 



p — 



dV 



T,{Nk} 



= ^ /ifcrife - Fnnc/V 



(4) 



(5) 



which is consistent with the definition of the grand thermodynamic potential, PV. We 
have to give these element ary details here since we find that some terms are omitted in the 
expressions for fj,i, e.g. bv lYakovlev fc Shalvbkov ( 1989t ). 

For example, let us start with the Coulomb contribution to Pmic- As usual we introduce 



with being the radius of the ion sphere: 



1/3 



47r 



''Zi 



(6) 



(7) 



If Qi denotes the Coulomb correction to the free energy and has the functional form 

g, = kBTUT,) (8) 



then we find the Coulomb contribution to the chemical potential as 



Coul 



knT 



(9) 



contrary to Eq.(17) of lMazurek et all (|l979f ) and to Eq.(97) of lYakovlev fc Shalvbko^ (|l989l ) 
who find /i^ = ksT f[j{Vi), but their expressions lead to zero correction to Pnuc in Eq. ([5]). On 
the other hand, our expression Eq. ([9]) gives the expected correction t o Pnuc after substitution 
into Eq. ([S]). Basically the formulae in lYakovlev fc Shalvbkov ( 19891 ) are correct, since one is 
free to separate Coulomb contributions, ascribing them either to electrons or ions, but one has 
to avoid confusion with this separation. In the book bv lHaensel et al.l (|2007l ) they mostly rely 
on free en ergy F instead of Mi - Actually, the expression which is equivalent to ^ appeared 



already in lDewitt et al. however, the argument on the ro le of the chemical p otentials 



and debates on approaches to them still persists. For exa mple, Brown et al. ( 2006| ) discuss 
the difference of the ir use of chemical potentials with that of Dewitt et al. 7x973 ). Essentially, 



(Il973l) for a 
tidl " 



Brown et al. ( 20061 ) found that the Coulomb correction found bvlDewitt et al. 
classical plasma is applicable to quantum plasma as well. (See also Brvdges fc Martinlll999l 
for a detailed review on modern per s pectiv e on Coulomb systems). 

19791 ) in the assumed form of the free energy but with 



We closely follow iMazurek et al. 
some modifications. The m ain results presented b elow use the Coulomb correction as in 
Eq.®, as was also taken bv iMazurek et al" ( 19791 ). In the end we check the effect of the 



additional term in Eq.®. Thus, we put 



W + Qp + Bp 



(10) 



(11) 
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Here and ^)^' are the Fermi-Dirac expressions for the free energy of non-relativistic 
nucleons and they are express ed through the Fer mi integrals of half- integer index. To calcu- 
late them we use the code bv iNadyozhinI ( 1974al |bh. (The superscript zero indicates that an 
expression is for non- interacting particles). Fo r the bulk inter a ction energy of free nucleons 
we take the zero temperature expre ssion from iMazurek et al.l (11979 ) where the y have used 
the expression of Bavm et al. ( 197l[ ) with minor corrections from iMackid (Il976l) : 



.(0) 



(12) 



with = 1.5 7r^(np + ?tn ) and = rir,/ (r i p+nn ). We have used the same values of parameters 
entering in Eq. ^ as iMazurek et all (l979l ): ko = 1.34 fm"\ Wo = 15.5 MeV; s = 27.1 
MeV; K = 268 MeV. 

For nuclei (i > 2) we assume: 



^^ = + Q^ + B, 

The binding energy in Eqs. HU]), ([II]) and is 

A, 



B, = M, 



An 



-Mr, 



(13) 



(14) 



with Mi being the mass of a nucleus number i, which depends on the choice of the reference 
nucleus m. We use ^^Fe as the reference nucleus in our calculations. In calculating the total 
energy we exclude electron mass irie from the electron energy, then we have 



B, = AM, 



Ar 



-AM„ 



(15) 



whe re AM, is the atomic mass excess, which accounts for the mass of electrons. To compare 
with lLattimer fc Swestvl (|l991r ) one has to take into account that they assume i?n = 0, where 

m = 2 = n in our notation, and Bi = AM^ — Ai/S.My^ 

Considering the Coulomb corrections we follow Lattimer et al.l ( 1985f ) and take for the 
total Coulomb energy of a nucleus 



e: 



Coul 



3 ZFe2 



[1 



3 rA }_('[±Y 



5 rA 

Here ta is the nuclear radius, va — 1.2 • 10^^'^A^^'^ cm. So the Coulomb correction is 



Q^ 



3 Zfe' 



5 rA 



3 TA 

2az 



-(-)'] 



(16) 



(17) 



In the leading term, which is what we need at high density, this is just —0.9Ti with F^ 
defined in Eq.(l6]) when the form ([8|) is used. With the modificatio n as in Eq.(|9t it wiU be 
— 1.2Fi - very close to the numbers in Eqs. (56) and (57) found by Dewitt et al.l ( 1973h by 
two different methods. 

The Boltzmann expression for the chemical potential satisfies the relation 



exp(/z, //cbT) = ^ 



h^No 



3/2 



^2nkBTA^ 

One has to introduce the corrections A^ in addition to nf: 

= Z/i^ + {A- Z)^i^ + A, , 

where 



(18) 



(19) 



A,; 



, Coul 



(20) 
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The corrections for degeneracy are: 



and for the Coufomb part: 

Af = ZQp - . 
For nuclear interactions we have: 



(21) 



(22) 



+ AW + (Xp + X„) 



rydW I I A 7\dW 



' dXp 
dXA 



dX„ 



X. 



dXjs 



Z'>1 



(23) 



Mazurek et al. ( 1979[ ) have argued that these expressions can often be simphfied. We 



take relation Eq. ([23]) in full which requires an additional loop of NSE iterations. 
We assume the following correction to the surface energy: 



w 



1/2 



l-(Xp+X„)- 



4/3 



(24) 



The expression actually used in the mass formula for the surface energy may be much more 
complicated than just Wgur^^''^, but since the term in braces tends to zero for vanishing p 
this correction does not influence the vacuum values of nuclear masses. 



2.2 Spin of ground state 

To find the spin of ground states of exotic and experimentally unknown nuclei we have used 
the simple shell model with the Woods-Saxon potential for defining the scheme of levels and 
nuclear systematics. It is well known from experiments and calculations using up-to-date 
versions of shell model with residual interaction that the spin of ground states of even-even 
nuclei is equal to zero. The spin of nuclear ground states of odd-^ nuclei is determined 
generally by the spin of the uncoupled nuclcon (exceptions: ^^Fg, ^'^Nan, ^^Mn25). That is 
why in the present paper we defined the ground state spin of odd- A nuclei by the spin of an 
uncoupled proton or neutron according to experimental data and shell model calculations 
with phenomenological coupling. 

The spin of ground state of odd-odd nuclei is defined by the spins of the uncoupled 
proton, ip, and neutron, j„, and could not be calculated precisely even on the base of the 
microscopic models of atomic nuclei. Th at is why we hav e used, in the odd-odd case, the 
simple phenomenological Nordheim rules (Nordheim 1950? ). 



In the case of opposite spins of proton and neutron (j„ = /„ ± 1/2, jp — Ip T 1/2) the 
spin of an odd-odd nucleus is equal to 

I ^ in- jp ■ (25) 

In the other case (j„ = Z„ ± 1/2, jp — lp± 1/2) it is rather difficult to find the right value 
of ground state spin because it can fall into the wide range: 

jn -jp<I< in + ip ■ (26) 
But it is well known from systematics that / « j„ + jp mainly. For simplicity we have chosen 

I = ip + in- 

The value of the ground state spin should not strongly influence the results of NSE calcula- 
tions because in our model the nuclei are mostly in the excited sta tes. We have done a calcula- 



tion u sing another approach of the ground state spin calculation (jEngelbrecht fc Engelbrecht 



19911 ) to check the influence. It is also based on the simple shell model level scheme, and 
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significantly differs from ours only in the case of odd-odd nuclei, for which the spin is defined 
as: / = 0p+jn + l)/2. 

We find that the calculated results depend weakly on the different approaches to ground 
state spin calculation only in the region with the significant amount of "heavy" nuclei, where 
the role of odd-odd nuclei rises. Since the most abundant nuclei are even-even and odd-A 
ones under NSE conditions, one can use any simple estimation of the ground state spin 
of nuclei under these conditions. For Tg < 10 and low p one should prefer our model of 
ground state spin evaluation, or an "exact" calculations of nuclear characteristics from first 
principles. 



As for the role of the excited states, we agree with the conclusions of lLiu et al.l ([20071) that 
the effect is not very large. They used a reliable method of nuclear partition function calcula- 
tion, based on the Fermi gas formula. Ho wever, the question dese r ves further investigation 
especi a lly with account of new data from Rauscher et al.l (|l997l ): iRauscher fc Thielemann 
( 2000l ): lRausche^(|2003l ). The used in fLiu et all (12007^ is related to an e nergy-dependent level 
density parameter with micros c opic correction to a nuclear mass model ( Rauscher et al.flOQ? : 
Rauscher fc Thielemann 2000l : Rauscher 2003 ). Our code may incorporate various options 
for the partition functions. 



2.3 Nuclear masses 



We adopt the table of nuclear masses fKT UY') bvlKoura et al.l (j2005f ) for our EOS calcula- 
tions. The prediction of nuclear masses bv lKoura et al.l (120 05') is based on the mass formula 
composed of macroscopic and microscopic terms (jKoura et al. 20001. which treat deforma- 
tion, shell and even-odd energies. It covers ^9000 species of nuclei (Z > 2 a nd N > 2) in the 
nuclear chart. We use the extended mass table containing '^20000 species (jKoural 120071 ) to 
cover the neutron- and proton-rich regions in the nuclear statistical equilibrium. The root- 
mean-square deviation from experimental masses is 666.7 keV. We e xamine the depen dence 
on the nuclea r mass by ado pting other tables of the nuclear masses bv lHilf et al.l(|l976h . The 
mass data by MoUer et al.l (fl995 ) will be implemented in the code. Those mass tables predict 
the experimental masses very well around the stability line in the nuclear chart, however, 
they provide different predictions in neutron-rich region away from the stability. It would be 
interesting to examine whether these differences appear in the composition of dense matter 
in supernova core. 



2.4 Shen equation of state 

We adopt the table of EOS by IShen et al. I (Il998al lbh to investigate the infiuence of mixture 



of nuclei on the equation of state. The Shen EOS is a set of thermodynamical quantities 
of dense matter under the wide range of density, proton fraction and temperature for su- 
pernovae. It has been widely used for numeric al simulations of co re-collapse supernovae 
( Sumivoshi et a l. 2005, 2006; Janka ct al. 2005; Burrows et al. lEooi) and other astrophys- 
ical phenomena (jRosswog fc Liebendorferl l2003i) . The e quation of state is calculated by 
the relativistic mean field theory (jSerot & Walccka"l986'), which is based on the relativis- 
tic Briickener Hartree Fock theory (Brockmann & Machlcidt 19 90|) and is constrained by 
the experimental data of neutron-rich nuclei (jSugahara fc Tokilll994( l. The nuclei in dense 
matter is described within the local density approximation assuming one species of nucleus 
surrounded by neutrons, protons and alpha particles in the Wigner-Seitz cell. The basic 
behaviour of the Shen EOS in s upern ova core as well as the comparison with the an other set 
of EOS bv lLattimer fc Swestvl (|l99ll ) can be found in lSumivoshi et al.l (|2004 l2005l) . 



2.5 Beta-equilibrium 

Although we fix the value of in the current calculations, we can impose the condition of 
beta equilibrium in the Saha approach. We add the formulation here for the applications to 
proto-neutron stars. 
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The condition of beta-equilibrium is governed by the relation 



Mp + Mo = Mn + 



(27) 



Therefore, one has to calculate chemical potentials of neutrinos and electrons in an additional 
loop of iterations. 

For cal c ulating lepton ch emical pot entials one can use simple relations derived by 
Nadvozhin ( 1974a 3) (see also Blinnikov .1987: Bli nnikov fc Rudzskv. .! 



dP 



T 



(28) 



Here rt is a fermion number in the unit volume, i.e., the fermion concentration, h is an anti- 
fermion concentration, n — n is a charge of the unit volume, for example, for the electrons 
and it is a lepton charge for neutrinos {g = 2 for electrons and g — 1 for neutrinos). It is 
possible to express the lepton chemical potentials /i through the radicals if n — n and T are 
given. In actual calculations it is faster to do this by Newton iterations. 



3 Results 



We report the numerical results in a set of selected conditions for density, p, electron fraction, 
Ya, and temperature, T. We show the properties of dense matter at typical conditions 
(p=10^°-10^^ g/cm^) during the core collapse, where the composition of nuclei and free 
protons are crucial to determine the electron capture. We remark that the present framework 
is generally applicable in the wide range of conditions, which are necessary for core-collapse 
supernova simulations. However, it is limited up to ~0.3 due to the coverage of nuclear 
mass models at neutron-rich region. The limited range of low T ~2 MeV and high p ~10^'^ 
g/cm"^ is due to the consideration of nuclear interaction and the convergence of iterations. 
We show, at first, the general behaviour o f equation of state in the current framework 



with the mass formula by iKoura et al.l (|2005f ) in Next, we co mpare the equation of 



state obtained by the current Saha-approach with that of Shen EOS (jShen et al. I ll998bh 



We examine the difference between the two nuclear mas s mod els 

()l979h 



with the results calculated in the formulation of Mazurek et al 



We also compare 
In i j3.3[ we show the 



composition of nuclei in the nuclear chart in the situations of supernova core. 



3.1 General behaviour 

We show, in Figs. [H [2l [3] and IH the global features of dense matter at Y^ — 0.316 for a wide 
range of density and temperature relevant to collapsing supernova cores. We have calculated 
also the cases with 1^ = 0.398 and Y^ = 0.473, which are not shown here. The values of Y^ 
are selected so as to match with the values in the original table of Shen EOS. We see very 
smooth variations of calculated quantities in the range of density and temperature shown in 
the plot. There are some wiggles in the plots of the average of mass number. A, and proton 
number, Z in Figs. [21 [3] due to the shell effects. One can recognise the appearance of nuclei 
at low density ~ 10"*^° g/cm^ at the lowest temperature. This is different from the Shen EOS 
as we will see in the next subsection. 



3.2 Comparisons among models 

We compare the current results with those in the Shen EOS to see the effect of multi- 
composition with respect to the one-species treatment, which has been used in most of the 
supernova EOS's. We examine also t he dependence o n the nuclear mass models by comparing 
the cases bv lKoura etHI (|2005h and lHilf et all (|l976h . We set here Y^ = 0.316 and T=2 MeV 
which roughly corresponds to the neutron-rich supernova core. This T value also corresponds 
to the lower T boundary in the figures shown above. 
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We show that the mass fractions of nuclei, free protons and alpha particles as a function 
of density in Figs. [5] and [6l As the density goes higher, we have found the nuclei appear 
earlier than the Shcn EOS does. This is because the nuclei other than alpha particles appear 
abundantly in the treatment of multi-composition. In the Shen EOS, on the other hand, 
alpha-particle appear more in this density region ~10^^ g/cm^ as seen in Fig. [SI Figure 
[7] elucidates clearly this difference further. From '-^lO^ g/cm'^, the Saha-treatment provides 
nuclei of average mass number A ~ 10. At densities higher than ~ 3x10^^ g/cm"^, the 
average mass number is smaller than that of representative nuclei in the Shen EOS. The 
one-species treatment may overestimate the mass number in this respect. We note that 
the average mass number in the Saha-treatment does not include the contribution of alpha 
particles, which are separately treated in the Shen EOS. 

The thermodynamical quantities, i.e. entropy and pressure of nuclear contribution, are 
shown in Fig. [51 The entropy of the Shcn EOS is in accord very well with the current result. 
The pressure of the Shen EOS has trough around 10^^ g/cm'^ because of the dec rease of 



nuclear part of free energy. As it has been already reported in I Shen et al.l (|l998a^ . this is 
due to the increasing binding energy of nuclei. We note that one should add lepton and 
radiation contributions in addition, therefore, the total pressure is always positive. 

In Fig. [21 we show the average mass number and proton number of nuclei for = 0.473 
and T=0.63 MeV, which is more closer to the condition at centre of the initial progenitors. 
We can see that mass number (proton number also) in the current treatment is larger at low 
density and smaller at high density. This is quite similar to the case of more neutron-rich 
and higher temperature, 1^ = 0.316 and T—2 MeV in Fig. [71 

We comme n t here on the difference from the calculated EOS using the treatment by 
Mazurek et al. ( 19791 ). The calculation done without the second line in Eq. (123]) . which 



express the dependence of binding energy of nucleus on the nucleon fractions, are shown by 
dash-dotted lines in Figs. [51 to [51 In general, the treatment without iterations gives similar 
results to the one with iterations to calculate consistently the modified binding energy due 
to the presence of nucleons at this condition. In Fig. [TOl we show the mass fraction of free 
pr otons as a funct i on of density for Fe = 0.398 and T=2 MeV. We have found the treatment 
bv lMazurek et al. ( 19791 ) provides similar results to those of the Shen EOS and our treatment 



gives lower proton fraction by large factors. This may influence the electron capture rate in 
collapsing supernova cores. 

We remark on the dependence of mass models. The calculations by two mass models 
are similar to each other and one can hardly see the difference. Only in the prediction of 
mass number and proton number, there are slight differences between the two lines. In 
this respect, the derivation of bulk quantities, such as entropy and pressure, can safely be 
calculated without the uncertainties of mass formulae. We note, however, that the detailed 
composition depends on the adopted mass formula as we will show below. 

3.3 Nuclear composition 

We explore the calculated composition using the mass formula bv iKoura et al. (2Qq3) in 



Figs. [TTl [121 and \W\ for three typical conditions of supernova core. We took the condition s 
from the numerical simulation of core-collapse from a 15Mq star bv lSumivoshi et al (20o3)- 
At p=10^^ g/cm^, where electron captures goes on during the gravitational collapse from 
the initial Fe core, nuclei up to A ~ 100 appear mostly with the peak abundance a,t Z 
34 and the magic number A^=50. By the time the density reaches p=10^^ g/cm^^, electron 
capture is about to cease and becomes smaller. Because of higher temperature due to 
the collapse, the distribution becomes wider, reaching A ~ 132 at the double magic nuclei, 
^■^^Sn. Neutron-richness (small Yc) makes the peak position at more neutron- rich but at 
iV=50. At high density of 10^^ g/cm'^ and high temperature of 3 MeV, which corresponds 
to the moment just before the core bounce, the distribution extends from nucleons to A ^ 
160 continuously. We see the effect of magic numbers for neutron number 50 and 82 and 
even-odd numbers, which exist in nuclear mass data. 
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In order to demonstrate the dmerence or com position d u e to t he mass models, we show 
the case of g/cm^ with the mass data bv iHilf et al. ( 1976[ ) in Fig. [TH Through the 



comparison with Fig. 1121 the distribution is narrower both in A'^ and A = N + Z directions. 
The peak position is quite close each other in the two models. We see similar differences 
also in the case of p=10^^ g/cm'^ and p=10^^ g/cm^, especially on the strength of magic 
numbers. Neutron-richness at high density causes larger differences among them. 

3.4 Difference in Coulomb corrections 

As we have discussed above, there is some controversy on the form of Coulomb corrections in 
the literature. We compare in Fig[Tn]the effect of changing corrections as in formula ([5]) and 
in on the proton fraction. The effect is visible only at highest density but still very small. 
A bit more pronounced effect is obvious in Fig |16l for the average mass of heavy nuclei again 
at highest density when the NSE model is near the border of applicabil ity. However, a more 
detailed study is needed here along the lines undertaken recently by iNadyozhin fc YudinI 
1 200^ . 



4 Summary and discussions 

We have studied the properties of hot and dense matter at the conditions relevant to su- 
pernova cores in the extended Saha approach. The multi-composition of nuclei is taken into 
account by solving the chemical equilibrium among nuclei, neutrons and protons. The con- 
tributions to chemical potentials due to the degeneracy of nucleons. Coulomb effects, nuclear 
effects, surface effects are included in the Saha formulation. The modification of binding 
energies of nuclei due to the existenc e of nucleons outside n ucle us is solved self-consistently . 
This is different from the treatment of lMazurek et all (|l979f ). Cf. iNadvozhin fc YudinI (|2005l ) 



were the difficult y of self-cons istent treatment is pointed out. We have adopted the recent 
mass formula by Koura et al.l 120 05). covering ^20000 nuclides, and compared the results 



with the ones bv lffilf et all (|l976[ ) 



We have calculated extensive grids of density and temperature for several choices of Yc- 
We found that the obtained quantities behave in general very smoothly. When we compared 
with the values in the Shen EOS, which is a popular set of EOS for supernova simulations, 
they coincide with each other closely for thermodynamical quantities. This may suggest 
that the one species treatment of nuclei adopted in most of modern supernova EOS's is a 
good approximation in this regard. However, the appearance of nuclei starts at lower density 
than in the Shen EOS and the mass fractions of alpha particles and nuclei are different in 
two codes. The average mass and proton numbers of nuclei are larger at low density and 
smaller at high density than those obtained from the Shen EOS. The former arises from 
the appearance of light nuclei around A ^ 10 in addition to alpha particles which are only 
included in Shen EOS. We explore also the composition in supernova core in the nuclear 
chart. We found that the distribution extends widely for high density and temperature and 
it depends on the choice of mass formula, which provides different shell effects. 

Our approach is rather simple and may be compared with more sophisticated modern 
tre atment of nuclear reac tions based on the statistical multi-fragmentation model (SMM): 



see 



Botvina fc Mishustin 

(|20M). Recently an exploration on the difference between the 
treatments under single and multi nuclear species approximations with phcnomenological 
mass formulae has been carried out to demonstrate the overestimation of mass number in 



the single nuclear species ISouza et al ] (jibos) 



The current finding assists understanding of the gravitational collapse of massive stars for 
supernova explosions. The difference of composition from the one species treatment in the 
Shen EOS suggests that the proper treatment in multi-composition of nuclei is important and 
may significantly affect the dynamics of core collapse. Since the electron capture on nuclei and 
free protons is crucial to determine the size of bounce core, i.e. the location of initial launch 
of shock wave, one should take into account the composition properly. Such efforts have been 
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made to evaluat e the total electron capture rates on nuclei by averaging the distributions 
(|Hix et alJ l2003l ) . One should, however, work carefully to predict the composition, since 
the mass and the nuclear effects may affect the detailed abundance. It is also necessary to 
calculate the equation of state, thermodynamical quantities and compositions in a consistent 
manner. We are aiming to apply the extended Saha treatment to provide EOS sets for 
supernova simulations in future. 
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Figure 1: The isocontours of the mass fraction of heavy nuclei, Xaz, (left) and free protons, 
Xp, (right) for Ye = 0.316 are shown in the plane of density and temperature. 
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Figure 2: The isocontours of the average mass number, A, of heavy nuclei, for = 0.316. 
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Figure 3: The isocontours of the average proton number, Z, of heavy nuclei, for = 0.316. 




Fi gure 4: The isocontours of the nuclear contribution to entropy, S^^c, for Ye = 0.316 in the 
plane of density and temperature. 
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Figure 5: The mass fraction of nuclei, Xa, (left) and free protons, Xp^ (right) as a function of 
density for 1^ = 0.316 and T = 2 M eV. The solid and das hed lines denote the NSE calculation 
with mass- formulae bv lKoural (|20o3) and lffilf et all (Il976l ^. respectively. The dash-dot lines are 
obtaine d using iKourai tmi) but with a simplified treatment of rel a,tion Eg. (|23|) as ha s been 
done by Mazurek et al. ( 19791 ). The dotted lines are the results from lShen et al.l ( 1998bl ). 




Figure 6: The mass fraction of alpha-particles, Xa-, as a function of density for 1^ = 0.316 and 
T = 2 MeV. The notation is the same as in Fig. [5j 
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Figure 7: The average mass number, A, (left) and average proton number, Z, (right) of nuclei 
as a function of density for = 0.316 and T = 2 MeV. The notation is the same as in Fig. [5j 




Figure 8: The entropy, S^uc (left), and pressure Pnuc (right), of nuclei as functions of density for 
Yf. = 0.316 and T = 2 MeV. The notation is the same as in Fig. [5j 
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Figure 9: The average mass number, A, (left) and average proton number, Z, (right) of nuclei 
as a function of density for = 0.473 and T = 0.63 MeV. The notation is the same as in Fig. [5l 




Figure 10: The mass fraction of free protons, Xp, as a function of density for = 0.398 and 
T = 2 MeV. The notation is the same as in Fig. [5j 
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p=10"gW, T=l MeV, Y^=0.398 
Mass: Koura et al. 
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Figure 11: The mass fraction of nuclei, Xa-, having proton number, Z, and neutron number, 
iV = A - Z in the nuclear chart for p = lO^^^/cm^, = 0.398 and T = 1 MeV. 
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p=10'^g/cm\ T=2 MeV, Y^=0.316 
Mass: Koura et al. 
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Figure 12: The mass fraction of nuclei, Xa, having proton number, Z, and neutron number, 
iV = ^ - Z in the nuclear chart for p = lO^^y/cm^, = 0.316 and T = 2 MeV. 
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p=10"g/cm\ T=3 MeV, Y^=0.316 
Mass: Koura et al. 
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Figure 13: The mass fraction of nuclei, Xa, having proton number, Z, and neutron number, 
N = A — Z in the nuclear chart for p = W^'^g/cm^, Y^, = 0.316 and T = 3 MeV. 
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Figure 16: The isocontours of the average mass number, A, of heavy nuclei, for = 0.473. Left 
paneh the Coulomb correction is taken as in Eq.®. Right panel: as in Eq.Q 
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